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ABSTRACT 

We have performed a cross-identification between OGLE-II data and single-epoch SIRIUS 
near-infrared (NIR) JHK survey data in the Large and Small Magellanic Clouds (LMC and 
SMC, respectively). After eliminating obvious spurious variables, variables with too few good 
data and variables that seems to have periods longer than the available baseline of the OGLE- 
II data, we determined the pulsation periods for 9,68 1 and 2,927 variables in the LMC and 
SMC, respectively. 

Based on these homogeneous data, we studied the pulsation properties and metallicity 
effects on period-ZT magnitude {PK) relations by comparing the variable stars in the Large 
and Small Magellanic Clouds. The sample analyzed here is much larger than the previous 
studies, and we found the following new features in the PK diagram: (1) vari able red giants 
in the SMC form parallel sequences on the PK plane, just like those found bv lWood ( I200Q,) 
in the LMC; (2) both of the sequences A and B of Wood ( 2000) have discont i nuitie s, and they 
occur at the /T-band luminosity of the TRGB; (3) the sequence B of "Wood" (2000) separates 
into three independent sequences B^ and C'; (4) comparison between the theoretical pulsation 
models ({Wood & Sebo 1996) and observational data suggests that the variable red giants on 
sequences C and newly discovered C' are pulsating in the fundamental and first overtone 
mode, respectively; (5) the theory can not explain the pulsation mode of sequences and 
B^, and they are unlikely to be the sequences for the first and second overtone pulsators, 
as was previously suggested; (6) the zero points of PK relations of Cepheids in the metal 
deficient SMC are fainter than those of LMC ones by ps 0.1 mag but those of SMC Miras 
are brighter than those of LMC ones by « 0.13 mag (adopting the distance modulus offset 
between the LMC and SMC to be 0.49 mag and assuming the slopes of the PK relations are 
the same in the two galaxies), which are probably due to metallicity effects. 

Key words: galaxies: Magellanic Clouds - stellar content - stars: AGE and post-AGB - 
variables - infrared: stars - surveys 



1 INTRODUCTION 

The Large and Small Magellanic Clouds offer us the opportunity 
to study the stellar evolution. They are close enough to resolve 
each stellar populations with relatively small telescopes and yet far 
enough to neglect their depth, so that we can reasonably consider 
that all stars in them are at the same distance from us. 

Many observations have been carried out toward the LMC 



and SMC, especially for the variable stars to study the relation- 
ship between the pulsation period and the luminosity {PL relation- 
ship). There is a long standing question related to the PL rela- 
tionship, whether it depends on metallicity or not. A controversy 
still exists between observations and theories (e.g., 'Alibert et alj 
1999 and r eferences ther e in for Cepheid variable s; Wood & Sebo| 
1996, Glass etal. 199^ iFeast. Whit elock & Menzies 2002 for 
long-period variables). Because the PL relationship is of fundamen- 
tal importance for the determination of the extragalactic distance 
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scale, a homogeneous survey of variable stars in different environ- 
ments is needed. 

Thanks to the gravitational microlensing search projects 
(OGLE, e.g Udalski. Kubiak & Szvmahski 1 997.; M ACHO, e.g., 
AlcocketaljEoqOl: EROS, e.g., lAfonso et al]fl999l : MOA, e.g., 
BondetalJl200lh . a large num ber of varia ble stars have been 
found in the Magellanic Clouds. IWood et alj fl999) analyzed the 
data from MACHO survey in the 0.25 square degree area of the 
LMC bar and showed that the variable red giants form paral- 
lel sequences in the period-magnitude plan e. ICioni et ^ i200ll) 
andjLebzelter. Schu ltheis & Melchio? ^2002*) confirmed the results 
of|Wood (2000) using data from EROS mic rolensing survey and 
AGA PEROS variable star catalog I Me lchior. Hughes & GuibertI 
|200Q^ . respectively, both in the 0.25 square degree area of the LMC 
bar. 

In this paper, we study the OGLE variables in the 3 square 
degrees along the LMC bar (see figure 1 of Ita et al. 2002) and in 
the 1 square degree area of the central part of the SMC (figure^. 
The OGLE variables have been cross-identified with the SIRIUS 
NIR survey data and their pulsation periods were determined by ap- 
plying the Phase Dispersion Minimization technique to the OGLE 
data. We classify the variable stars on the PK plane, and discuss 
the properties of each group. The environmental effect on PK re- 
lations is discussed by comparing variable stars in the Large and 
Small Magellanic Clouds. 
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Figure 1. A schematic map of the monitoring region in the SMC. Each mesh 
has a side of 20 arcmin and consists of 9 sub-meshes, which correspond to 
the field of view of the SIRIUS. The background stars are drawn from the 
GSC2.2 catalog. Stars brighter than 16th magnitude at Fgsc located in the 
total area of 25 square degree centered on (a2000, 82000) = (o'^ 55™ 00.0^, 
-70° 30' 00.0") are plotted. 



2 DATA 

2.1 IRSF/SIRIUS 

We are now conducting a NIR JHK monitoring survey toward the 
Magellanic Clouds using the InfraRed Survey Facility (IRSF) at 
Sutherland, South African Astronomical Observatory. The IRSF 
consists of a specially constructed 1.4m alt-azimuth telescope to 
which is attached a large-format 3-channel infrared camera, SIR- 
IUS (Simultaneous three-colour InfraRed Imager for Unbiased Sur- 
veys). The SIRIUS can observe the sky in J H and Ks bands simul- 
taneously and has a field of view of about 7.7 arcmin square with 
a scale of 0.45 arcsec/p ixel. Details of the ins trument are found in 
iNaeashima et alJJl999l) and lNagavama et alj f2002). 

We have been monitoring a total area of 3 square degrees 
along the LMC bar since December 2000, and a total area of 1 
square degree around the SMC center since July 2001. Figure Q 
shows the monitoring area in the SMC . Details and so me first re- 
sults of the LMC survey are found in llta et alj i2002h (hereafter 
referred to as Paper I). In the present paper we deal with the single- 
epoch data from the monitoring survey. 



2.1.1 Photometry 

All data were reduced in the same way using the SIRIUS pipeline 
(Nakajima, private communication) and a single image comprises 
100 dithered 5 sec exposures. Photometry was performed on 
the dithered images with DoP HOT in the fixed position mode 
iSchechter. Mateo & Saha"l993\ The magnitudes are referred to 
the Las Campanas Ovservatory (L CO) system ones bas ed on obser- 
vations of a few dozen stars from lPersson et all il998h . In the sys- 
tem conversion, colour terms have been considered by observing 
red stars in table 3 of Persson et al. The photometry has a signal-to- 
noise ratio (S/N) of about 20 at 16.9, 16.8 and 16.0 mag in J, H and 



K, respectively. Typical photometric errors are tabulated in table 1 
of Paper I. 

2.1.2 Source selection 
LMC: 

In the three square degree survey area, 1,593,248 sources were de- 
tected by DoPHOT. Because of the scan overlaps, there should be 
many multiple entries of a single source in the original sample. 
We eliminated such multiple entries based on (I) spatial proximity 
(|Ar| 0.5") and (2) difference in photometry (\Amj [fg„jf[\ 0.1 
mag). It leaves 1,305,983 sources. Further, requiring detections at 
least two of the three wave bands leaves 789,929 sources. 

In Paper I, we used the allowed discrepancy \Ar\ of 1.5", in- 
stead of the present value of 0.5". This was because the celestial 
coordinates were calculated by referencing the USNO catalog in- 
stead of the Guide Star Catalog II (GSC2.2, see section 2.1.3.). We 
can find substantial number of positional reference stars by using 
the GSC2.2 catalog, and this greatly improved the accuracy of the 
positioning. 
SMC: 

In the one square degree survey area, 257,637 sources were de- 
tected by DoPHOT. Multiple entries were eliminated in the same 
manner as we did in the LMC. It leaves 210,566 sources. Further, 
requiring detections at least two of the three wave bands leaves 
107,376 sources. 

2.1.3 SIRIUS data 

In table □ we summarize the number of sources detected in each 
wave band. Hereafter, we refer these 789,929 and 107,376 sources 
as the SIRIUS data. Although many stars were detected only in 
two wave bands, these stars hardly affect the results of this paper. 
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Table 1. Breakdown list of the SIRIUS data 







Number 


Wavebands 




LMC 


SMC 


/ + // + 


K 


596,913 


85,491 


J + H 




75,812 


15,173 


H + 


K 


114,827 


6,693 


/ + 


K 


2,377 


19 


Total 




789,929 


107,376 



Also, although SIRIUS JHK photometries are single-epoch ones, 
it should not affect the results so much, because the pulsation am- 
plitudes of variable stars are usually small in the infrared (except 
for Mira variables, but still the results should be statistically fair 
because of large number of available data). 

The SIRIUS data were corrected for the int erstellar absorp- 
tion based on the relations in Koornneej il982h . assuming R = 
-i.l. We adopted (Ay, Ah, Ak) = (0.129, 0.084, 0.040) and 
(0.082, 0.053,0.025), corresponding to the total mean reddenings 
ofEB-v = 0.137 and 0.087 as derived by Udalski et al. ( 1999a) for 
the ogle's observing fields in the LMC and SMC, respectively. 



2.1.4 Astrometry 

Celestial coordinates of the sources detected by the SIRIUS were 
systematically calculated in the International Celestial Reference 
Syste m (ICRS) by referencing the Guide Star Catalog II iSTScj 
I2OOII) . The astrometric precision of the SIRIUS data is expected to 
be the same as that of GSC 2.2 catalog, i.e., better than 0.5" in most 
cases (Tanabe et al. in preparation). 



2.2 OGLE 

Thanks to the microlensing surveys (OGLE, MACHO, EROS, 
MOA), a large number of variable stars were found as the natural 
by-products. In particular, the Optical Gravitational Lensing Ex- 
periment (OGLE) project has been monitoring the central parts of 
the Magellanic Cloud s in the BVI wave b ands and the / band time- 
series data (OGLE-II;'Udalski et al."l997') is now available over the 
Internet (OGLE homepage; http://sirius.astrouw.edu.pl/""ogle/l. 



2.2.1 Cross-identification of the SIRIUS sources and the OGLE 
variables 

The SIRIUS data and all the stars found by the OGLE survey in 
the Magellanic Clouds have been cross-identified using a simple 
positional correlation. In the very strict sense, the OGLE survey 
does not fully cover the SIRIUS survey area, and some very small 
areas are not included in the OGLE. However, it never affect any 
conclusions in this paper. 
LMC: 

There are 49,008 OGLE variables in the OGLE/SIRIUS overlapped 
3 square degree area in the LMC. The identification was made in 
two steps. First a search radius of 3" was used, finding 35,726 
matches. Plotting the offsets between the SIRIUS coordinates and 
the OGLE coordinates revealed that there are small median astro- 
metric shifts between them; the median differences are -0.80" in 
right ascension (R.A.) and 0.10" in declination (DEC). Second, the 
OGLE coordinates were corrected for these offsets, and a search ra- 
dius of 3" was used again. We finally found 35,783 matches with an 
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Figure 2. Histogram of the positional difference between the SIRIUS co- 
ordinates and the OGLE coordinates in 0.05" bins for the cross-identified 
stars in the LMC (upper panel) and the SMC (lower panel). 



rms error in the difference between SIRIUS and OGLE coordinates 

0.68" in R.A. and 0.63" in DEC, respectively. 

SMC: 

There are 7,345 OGLE variables in the OGLE/SIRIUS overlapped 
1 square degree area in the SMC. The identification was made in 
the same way as we did in the LMC. First a search radius of 3" 
was used, finding 6,118 matches. The offsets between the SIRIUS 
coordinates and the OGLE coordinates are 0.02" in R.A. and -0.40" 
in DEC. After correcting the OGLE coordinates for these median 
offsets, a search radius of 3" was used again. We finally found 6, 103 
matches with an rms error in the difference between SIRIUS and 
OGLE coordinates 0.41" in R.A. and 0.38" in DEC, respectively. 

We show the histogram of positional differences of the finally 
identified stars in Figure|2| It is clear that the positional differences 
are smaller than 1" for most of the identified stars (29,228 out of 
35,783 stars or 81.7% in the LMC, and 5,745 out of 6,103 stars or 
94.1% in the SMC). 



2.2.2 Period determination 

We applied the Phase Dispersion Minimization (PDM) technique 
to the OGLE data to determine the light variation periods of 
the cross-identified stars. This technique is described in detail by 
IStellingwer3h978h and here we show only the outline of the tech- 
nique we used. 

We represent the observational dataset as (x,-,/,), wherexis the 
magnitude, / is the observation time and ; denotes the rth observa- 
tion. We assume that there are N observations in all (/ = 1,^). Let 
V be the variance of the magnitude of the original dataset, given by 



N 

/=l 



(1) 



where x is the mean; First, the x; and f; are folded with a 

trial period and the full phase space [0, 1) is divided into an appro- 
priate (user-specified) number of bins (there are Nb bins in all, i.e., 
i = 1,Nb). Second, the variance vj of magnitude of each bin and 
the mean variance v of them are calculated, namely. 



Nb 
7=1 



(2) 
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Table 2. Used resolution in period search 



trial period [day] 

startCJS) end(<) resolution [day] 



0.1 


1.0 


0.00001 


1.0 


5.0 


0.0001 


5.0 


10.0 


0.001 


10.0 


100.0 


0.01 


100.0 


1000.0 


0.1 



Third, the statistic parameter 9, which is defined as 

e = v/V, (3) 

is calculated for the trial period. These three processes are then 
repeated for the next trial period (i.e., P^^i^il new ~ Arial old + 
resolution; see below). The 6 is a measure of the regularity of the 
light variation or the scatter about the mean light curve, being hope- 
fully near for the very regular light variation with true period and 
near 1 for the irregular light variation without any distinctive pe- 
riodicity or for the incorrect period. Hereafter, we call stars with 
small 9 0.55) as "regularly pulsating variables", and the others 
(> 0.55) as "less regularly pulsating variables". 

We searched the minimum of the 9 throughout the trial peri- 
ods that start from 0. 1 and end in 1,000 days. The trial periods were 
incremented by variable step values (= resolution), which are sum- 
marized in tableU For all of the 35,783 (LMC) and 6,103 (SMC) 
matched stars, their original light curve, 9 spectra and the folded 
light curve are carefully eye-inspected, and the obvious spurious 
data points were eliminated. We assigned one pulsation period to 
each variable star, and the multi-periodic stars (e.g.. Bedding et al. 
Il998 f) were not analyzed in this paper. It should be noted that the 
total available baseline of the OGLE-II data is about 1,200 days and 
we did not try to determine the light variation period for the vari- 
ables that seem to have the period longer than 1,000 days. Also, 
we did not try to determine the period for the variables that have 
too few good data points and that seem to hardly show any pe- 
riodicity. The OGLE-II survey data has yearly blank phases that 
could be a problem for variable stars with pulsation period of about 
a year. However, experimentations showed that these yearly gaps 
hardly affect the period determination, and it is the virtue of the 
PDM technique. In very rare cases, the yearly gaps may unfortu- 
nately coincide with the timing of the peak luminosities. They can 
affect the determination of pulsation amplitude, but still, such cases 
should be very rare. Finally, we determined the light variation pe- 
riods for 9,681 and 2,927 OGLE variables in the LMC and SMC, 
respectively. 



3 DISCUSSION 

3.1 K magnitude distribution 

Many low-amplitude variables were newly discovered by the afore- 
mentioned optical surveys. The vast majority are at the luminosities 
below the TRGB. lAlves et al] 111998^ suggested that they are AGB 
stars in the early evolutionary phase prior to entering the thermal- 
pulsing phase of the AGB. In Paper I, we found that the K mag- 
nitude distribution of the variable red giants in the LMC has two 
peaks; one is well above the TRGB, consisting genuine AGB vari- 
ables and the other is just around the TRGB. 

In figure|3| we show the K magnitude distribution N{K) (dK = 
0.05 mag) of the SIRIUS data. The discontinuity in the N{K) is 
clearly seen around K ^ 12.1 (upper panel) and ~ 12.7 (lower 
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Figure 3. The K magnitude distribution differentiated by 0.05 mag bins of 
the SlRlUS data in the LMC (upper panel) and the SMC (lower panel). 
Closeup around the TRGB is shown in the inset. 

panel), which corresponds to the TRGB of the LMC and SMC, 
respectively. Figure|4]is the same diagram as the figure|3] but the 
9,681 and 2,927 OGLE/SIRIUS variables in the LMC and SMC are 
shown. It is clear that there is a "pile-up" around the A^-band lumi- 
nosity of the TRGB in both two galaxies. In Paper I, we concluded 
that a significant fraction of the constituent variable stars of the sec- 
ond peak could be on the RGB. Of course there may be faint AGB 
variables, but there is no reason to assume, according to models, 
that they accumulate at the TRGB. 

3.2 Period-A" magnitude diagram 

IWood i200Cl) found that radially pulsating red giants in the LMC 
form "three" distinctive parallel sequences in the PK plane. He 
characterized the complex variability of the red giants by the 
pulsation modes, whic h are an im p ortant probe of the internal 
structure of the stars. IWood etaP (,1999') suggested that Mira 
variables are fundamental mode pulsators, while so-called semi- 
regular variables can be pulsating in the first, second or even 
higher overtone mode, by comparison of observations with the- 
oretical models. However, a controversy still exists on the pul- 
sation mode of Mira variables. Diameter measurements of Mira 
vari ables in the Milky Way indicates first overtone puls ation (e.g., 
Hani ff. Scholz & TuthilJl995lWhiteIock & Feasj2o6ol) . while ob- 
servations of radial velocity variation favor fundamental mode pul- 
sation for Mira variables (e.g., Bessell, Scholz & Wood 1996). 

We show the relationship between the determined pulsation 
periods and the K magnitudes of the cross-identified variables in 
figure |5| We found that the variable red giants in the SM C form 
parallel sequences on the PK plane, just like those found by I Wood 
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Figure 4. The K magnitude distribution differentiated by 0. 1 mag bins of 
the OGLE/SIRIUS 9,681 and 2,927 variables in the LMC (upper panel) and 
the SMC (lower panel), respectively. 
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Figure 5. Period-S" magnitude diagram for OGLE variables in the LMC 
(upper panel) and the SMC (lower panel). 



teOO(f) in the LMC. In the LMC PK diagram, bunch of variable 
stars exist in the faint magnitude {K < 15.0) and very short pe- 
riod (logP < —0.1) region. Considering their location on the PK 
plane, they are likely to be RR Lyrae variables. Compared to the 
LMC PK diagram, the SMC one is a bit scattered. This is probably 
due to the difference in the intrinsic depth of the two galaxies, be- 
cause the LMC has a nearly face on orientation ( Westerlund 1997), 
and Groenewegen ( 2000) obtained the total front-to-back depths of 
0.043 mag for the LMC and 0. 1 1 mag for the SMC, respectively. 

In figure|6| we again show the PK relationship of variable stars 
in the Magellanic Clouds. In this diagram, we divided variable stars 
into four types according to their J — K colours and as is shown on 
the upper-left comer of the figure. The labels of the sequences are 
named in analogy to the ones found in the LMC bv lWood i200(t) 
(except C', F and G). While Cioni et al. ( 2001) found no objects on 
sequence A from their EROS/DENIS variables in the LMC, this se- 
quence is clearly seen in both galaxies. This is essentially the same 
diagram as the one obtained bv lWood J2000h . However, the sample 
used here is much larger than those of previous studies, and we can 
see new features in the figure: (1) the sequences A and B of Wood' 
^2000 ) are found to be composed of the upper and lower sequences 
A+,5+, and A",S"; (2) the sequence B of iWood (.2000) separates 
into three independent sequences and C'; (3) radially pulsating 
red giants form actually "four" parallel sequences {A,B,C' and C). 



3.2.1 Less regularly pulsating variables 

Most of the less regularly pulsating variables reside in the se- 
quences A and B, and show small pulsation amplitudes (A/ < 0.2 
mag, but some of the B stars have AI of about 0.4 mag). Al- 



though it is not mentioned in IWood<200(i) . his sequences A and B 
are densely populated at fainter magnitudes, and a discernible gap 
around the TRGB can be seen in his diagram (see his figure 1). This 
gap also presents in our diagram around the AT-band luminosity of 
the TRGB {K ~ 12.1 in the LMC and ~ 12.7 in the SMC), which 
we roughly estimated in section 3.1. Also, in the LMC, a horizontal 
misalignment of the sequences apparently occurs at the TRGB. One 
possible explanation for the discontinuity is that neither sequences 
A nor B are single sequences, but each of them is composed of two 
independent sequences. So we subdivided the sequences A and B 
into A+,S+ and A^,S^, where the + and — correspond to the lu- 
minosity, being brighter than K= 12.1 (LMC) or 12.7 (SMC) and 
fainter than it, respectively. 

Since the sequences A+ and 6+ exceed the A'-band luminosity 
of the TRGB, they consist of intermediate-age population and/or 
metal-rich old population AGB stars. However, the interpretation 
of the sequences A^ and B^ is difficult. If the metal-poor and old 
AGB stars pulsate in the same mode as the more massive stars on 
the sequences A+ or S+, they should be located on the sequences 
A^ or B^ because they do not exceed the TRGB luminosity. The 
question is whether the sequences A^ and B^ contain RGB pul- 
sators or not. If AGB stars alone make up the A^ and B^ as well 
as A+ and fi+ sequences, the small gap and the horizontal shift be- 
tween them needs to be explained. If A^ and B stars, at least some 
of them, are pulsating on the RGB, then the discontinuity between 
the — and + sequences is understood in terms of different evo- 
lutionary stages. If we could apply the same period-temperature 
relation as that for Mira variables iAlvarez & Mermessieilll997h 
to the A^ and A+ stars, the period shift of SlogP = —0.083 be- 
tween the two sequences at the TRGB (K = 12.1, see table|3j cor- 
responds to the temperature difference of SlogTgg^ ~ 0.015. This 
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Figure 6. Period-X' magnitude diagram for OGLE variables in tlie LMC 
(upper panel) and the SMC (lower panel). Variable stars are classified into 
four types according to their J — K colours and 9. The labels of the se- 
quences are named in analogous to the sequences found bv iWooA klOOO.) in 
the LMC (except C', F and G, see text). 

is in good agreement with the temperature difference between the 
AGB and RGB stars with the TRGB luminosity, predicted by the 
stellar evo lutionary model for stars with 1~2 M,;, in the LMC 
iCastellani et al. 2003) . Though we prefer the RGB interpretation, 
it should be noted that a definitive answer has not been obtained. 

3.2.2 Regularly pulsating variables 

Figure Q shows the comparison between the observational period- 
K magnitude sequences of pulsati ng red giants in the L MC and 
the theoretical models calculated bv lWood & S ebo ( 1996). The lo- 
cation of sequenc e C is consistent with the LMC Mira sequence 
jpeast etalJl989l) . After classifying stars by their regularity of light 
variation, we discovered a new sequence C'. Most of the stars on 
sequence C' are regularly pulsating just like those on C, but their 
amplitudes are smaller than those of the C stars. We note, however, 
that the amplitudes of stars on sequence C' are larger than those 
of stars on B. The observed periods, luminosities and period ratio 
of sequences C and C' agree very well with those of the theoreti- 
cal f undamen tal (Pq) and first overtone mode (Pi ) pulsation models 
& Sebo 1996). Also, the theory predicts that overtone pul- 
sators have smaller amplitudes than fundamental ones. Therefore, 
all of the above observational facts and the agreements with the 
theory could be naturally explained if we consider that the stars of 
sequences C and C' are Mira variables pulsating in the fundamental 
and first overtone mode, respectively. However, the pulsation the- 
ory cannot reproduce the period ratio of the observed sequences A 
and B. To explain their pulsation modes, a new theoretical model is 
needed. 
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Figure 7. Period-^T mag nitude relations of vari able red giants compared 
with theoretical model of Wood & SebcHl996li . Their models for stars of 
masses of 0.8, 1 .0 and 1 .5 Mv, pulsating in the fundamental (Pq, they forced 
theoretical fundam ental mode period to fit the observed Mira relation of 
iFeast et alj <1989ft (dot-dashed line)) and first (Pi), second (P2) and third 
(P3) overtone modes are shown. Colours of the points are as in figure l6l 



It is interesting that, with only a few exceptions, all of the red 
giants with red colours {J — K > 1 .4) reside in the brighter magni- 
tudes of sequences C and C' . According to lNikolaev & Weinberg! 
1200Q) . stars with J — K > 1.4 in the LMC are primarily carbon- 
rich AGB stars, being considered as the descendants of oxygen- 
rich AGB stars. We subdivided stars on sequences C and C' into 
two subgroups on the basis of their J — K colour, being redder than 
1.4 and the others. Table|3|(see section 3.3) infers that those two 
populations may follow the different PK relations, and the carbon- 
rich (J — K > 1.4) stars could define a line of somewhat shallower 
slope with brighter zero point than that of the oxygen-rich stars, 
[peast et al. 1 1989) obtained PK relations for carbon- and oxygen- 
rich stars in the LMC. The calculated slopes and zero-points are 
in good agreement with their results. In figure|8| we plotted colour- 
magnitude and colour-colour diagram of probable carbon-rich stars 
on sequence C and C'. Interestingly, all of the carbon-rich stars with 
J — K colour redder than 2.0 are on the sequence C, meaning that all 
of them are pulsating in the fundamental mode. However, it should 
be emphasized that the obscured oxygen-rich stars such as OH/IR 
stars could have redder colours. So, it could be dangerous to as- 
sume that all of the stars with redder colours are carbon-rich stars. 
Infrared spectroscopic work would provide the definite conclusion 
on these issues. 

The sequence D is thought to be populated by the longer pe- 
riod of a binary system 1 Wood 2000). However, the explanation of 
this sequence is not clear yet. The periods of the group D stars are 
very long an d more long- t erm o bservations are needed for the clear 
explanation. IWood etal] ^99^ suggested that stars on the loose 
sequence E are contact and semi-detached binaries. Interestingly, 
their distribution seems to extend to the TRGB. Stars on sequences 
F and G are very regularly pulsating and have periods ranging 
from less than 1 day to more than 30 days, suggesting that they 
are Cepheid variables pulsating in fundamental and first overtone 
mode, respectively. 
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Figure 8. Colour-magnitude (upper panel) and colour-colour (lower panel) 
diagram of the carbon-rich stars on sequence C (crosses) and C' (dots). Two 
stars that are on sequence C are excluded from the colour-colour diagram, 
because they are outside the range. 

3.3 Period-A' magnitude relations 

In order to determine the PK relations, we defined boundary lines 
as in figure |9| that divide the variable stars into nine (LMC) and 
eight (SMC) prominent groups (sequence E is excluded). For 
groups C, C' and D the slopes of the boundary line s were chosen 
to be -3.85, which is obtained bv lHughes & Wooj (1990) for the 
slope of the PK relation for the blend populations of carbon- and 
oxygen-rich Miras in the LMC. The slope of the right side of group 
S+ was also chosen to be -3.85. For groups A^B^, and FG, 

the slopes of the boundary lines were chosen to be -3.35, -3.45 and 
-3.36, respectively, after the rough estimation of the slopes of each 
sequence. The width of the each box was chosen so that each of 
them contains the most probable body of a sequence and it mini- 
mizes the miss-classifications as much as possible. 

Because the scattering is rather large due probably to the 
SMC's intrinsic front-to-back depth, we did not tried to separate 
C' stars from fi+ ones in the SMC. Also, we did not calculate the 
PK relations of variable red giants in the SMC. The dotted lines in 
figure|9|are the least square fits of linear relation to each group, and 
table|3|summarizes the calculated PK relations. 

3.3.1 Metallicity effects on the period-K magnitude relations 

iGascoiend i\9l4\ empirically obtained the period-visual magni- 
tude relation for Cepheid variables and found that Cepheids in 




Log P [days] 

Figure 9. Period-A' magnitude diagram for OGLE variables in the LMC 
(upper panel) and the SMC (lower panel) with classification boxes. The 
dotted lines ai'e the least-square fits of linear relation to stars in each box. 

the metal deficient environment would be fainter than those with 
the same periods in the metal rich one. Based on this result, the 
Cepheids in the metal-deficient SMC should be fainter than those 
in the LMC if compared with the same period. On the other hand, 
|Wood 1 1990) suggested that if there was a factor of 2 difference 
in metal abundance between LMC and SMC long period variables 
(e.g., Miras), the LMC stars were fainter by 0.13 magnitudes in K 
than those in the SMC with the same period. Therefore, Cepheids 
and long period variables are expected to express completely dif- 
ferent reactions to the metal abundance. 

One possible way to test this anticipation is to compare the 
PK relations in the LMC and SMC, and it is done in figure [Tol 
The distances to the Magellanic Clouds are still uncertain, and 
we can only say with confidence that the SMC lies some 0.4-0.6 
mag beyond the LMC iColell99^ . We adopted the dis tance mod- 
uh to the LMC and SMC as 18.40 ^Nelso n et alfcood ) and 18.89 
( Harries. Hilditch & Howarth 2003), respectively. These distance 
moduli were obtained by using the eclipsing binary systems in the 
Magellanic Clouds, and the technique, while not entirely geometri- 
cal (reddening estimates are required), is free from the metallicity- 
induced zero point uncertainties. So, we shifted the K magnitudes 
of the SMC stars by —0.49 magnitude to account for the different 
distance moduli of the SMC and LMC. The red dots represent the 
LMC stars and the green dots correspond to the SMC stars. 

The figure shows, just as expected, that Cepheids in the SMC 
are fainter, but red giants in the SMC are brighter than those in the 
LMC if compared with the same period. Also, it is likely that the 
slopes of the PK relations are almost the same in the two galax- 
ies and only the zero points seem to differ. To estimate the dif- 
ferences quantitatively, we fixed the slopes of the PK relations of 
SMC Cepheids and Miras to the corresponding LMC ones, and cal- 
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Table 3. Period-X' magnitude relations for variable stars in the Magellanic 
Clouds of the for m K = ax logP + b . The K magnitudes are referred to 
the LCO system <Perssonetaljil99 8Y The a is the standard deviation in 
K, and N is the number of stars being used in the least square fitting after 
applying 3.0 a clipping. Stars on sequences C and C' in the LMC were 
subdivided into two subgroups according to their J — K colour, and their 
period-A' magnitude relations were also calculated. 



Group 


a 


h 


a 


N 


J-K 


LMC 


A- 


-3.284±0.()47 


17.060±0.065 


0.114 


1142 


All 


A+ 


-3.289±0.047 


16.793±0.077 


0.125 


510 


All 


B- 


-2.931±0.057 


17.125±0.091 


0.100 


584 


All 


B+ 


-3.356±0.052 


17.634±0.099 


0.160 


502 


All 


C 


-3.768±0.023 


18.885±0.046 


0.110 


693 


All 


c 


-3.682±0.109 


18.720±0.245 


0.122 


135 


>1.4 


c 


-3.873±0.031 


19.081±0.060 


0.105 


558 


^1.4 


c 


-3.520±0.034 


19.543±0.082 


0.198 


975 


All 


c 


-3.369±0.099 


19.165±0.250 


0.200 


447 


>1.4 


c 


-3.589±0.055 


19.698±0.126 


0.196 


528 


^1.4 


D 


-3.635±0.078 


21.718±0.207 


0.198 


472 


All 


F 


-3.188±0.019 


16.051±0.013 


0.095 


540 


All 


G 


-3.372±0.041 


15.574±0.016 


0.091 


317 


All 


SMC 


F 


-3.350±0.020 


16.722±0.010 


0.139 


606 


All 


G 


-3.280±0.()35 


16.133±0.009 


0.147 


403 


All 
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Figure 10. Feriod-A" magnitude diagram of variable stars in the LMC (red 
dots) and SMC (green dots). The K magnitudes of SMC stars are shifted 
by —0.49 magnitude (see text) to account for the difference in the distance 
moduli between the LMC and the SMC. 



culated the zero points. The fixed-slope solutions yielded the zero 
points of 19.409±0.018, 16.170±0.006 and 15.659±0.007 mag for 
sequences C, F and G in the SMC, respectively (after shifting their 
K magnitudes by —0.49 mag). Comparing these zero points with 
those of the LMC Cepheids and Miras, Cepheids in the SMC are 
fainter by ~ 0. 1 mag and Miras in the SMC are brighter by ~ 0. 13 
mag than the LMC ones with the same period. 



4 SUMMARY 

We cross-identified the OGLE-II survey data and the SIRIUS NIR 
data in the Magellanic Clouds. Based on the large and homoge- 



neous data , we studied the pulsation properties of variable stars 
in the Magellanic Clouds. We discovered a new sequence on the 
period-A^ magnitude plane, and identified it as the one for Mira 
variables pulsating in the first overtone mode. We also found that 
the variable red giants in the SMC form parallel se quences in the 
period-A" magnitude plane, just like those found bv lWooj i200Ci) 
in the LMC. The period- A" magnitude relations were derived, es- 
pecially for both fundamental mode and overtone Cepheids in the 
Magellanic Clouds. These results are obtained with the greatest 
numbers of available stars and most deep photometric data in the 
infrared up to now. Metallicity effects on period-A^ magnitude re- 
lations are studied by comparing variable stars in the Magellanic 
Clouds. The result suggested that Cepheids in the metal-deficient 
SMC are fainter than the LMC Cepheids, while Miras in the SMC 
are brighter than the LMC ones if compared with the same period. 
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